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Abstract 

A general calculation of Casimir energies — in an arbitrary number of dimensions — for 

qj I massless quantized fields in spherically symmetric cavities is carried out. All the most com- 

,S^ • mon situations, including scalar and spinor fields, the electromagnetic field, and various 

boundary conditions are treated with the uppermost accuracy. The final results are given as 

- analytical, closed expressions in terms of Barnes zeta functions. A direct numerical evalua- 

rxj I tion of the formulas is then performed, which yields highly accurate numbers of, in principle, 

S I arbitrarily good precision. 
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1 Introduction 



0| and, also, the influence of different boundary 



Calculations of Casimir energies in spherically symmetric situations have attracted the interest 
of physicists for well over thirty years now. Since the calculation of Boyer |l|], who computed the 
Casimir energy for a conducting spherical shell and found, to his surprise, a repulsive force, many 
different situations in the spherically symmetric context have been considered. For example, 
dielectrics were included |||] (for the case of plane, parallel surfaces see [Q|) and used later 
for possible explanations of sonoluminescence |01~ 01. Moreover, enormous interest has been 
attracted by the MIT bag model in QCD | 
conditions has been considered in detail 

Different methods have been used for dealing with the Casimir effect. Whereas in the earlier 
times the Green function formalism was preferred, in recent years different approaches — which 
make use of contour integral representations of the involved spectral sums — are commonplace. 
Although the idea for this method, in the specific context of Casimir energies, goes back to the 
early days of the subject ||2^, a systematical, effective and simple application of this approach 
in various contexts has only recently been achieved [^, 22, 23, |2^, ^, 27, 28|. 




The spectral sum which actually appears in the calculation depends on the regularisation 
used and may include a cutoff function, to dampen high frequency contributions [28| or, as in 

, complex powers of the eigenvalues pi 



22, 26, 



As a 



the zeta regularisation technique | 

result, the details of the computation may differ slightly, e.g., in the specific integration contour 

chosen, but all of them share the elegance of this method. 

In recent contributions we have further developed the zeta function technique, in combination 
with several contour integral representations. Given the deep connections among zeta functions, 
heat kernels and functional determinants |31]- [p^], one advantage of the method is that it can be 



applied, alternatively, to the calculation of heat kernel coefficients [26|, functional determinants 
M, M (see also 



55|, ^gl), as well as Casimir energies |2C]- [22|. This clearly shows that zeta 
functions serve as a unified framework in different areas of interest. 

Here we want to pursue this idea, by using the zeta function framework in a precise analysis 
of the Casimir energy as a function of the dimension of space. Previously it had been shown, that 
arbitrary space dimension can be treated elegantly by making use of Barnes zeta functions, where 
the dimension can be considered as a parameter |27, ^] . This has been applied to the calculation 
of heat-kernel coefficients and determinants and it will be here used to study the Casimir energy. 
Apart from dealing with arbitrary dimensions, we will introduce scalars, spinor fields and the 
electromagnetic field in a unified way, including the effects different sets of boundary conditions 
have on them. In spirit, our analysis is to be compared with the one of Ambj0rn and Wolfram 



in Ref. |38], with the difference that the role of the Epstein zeta function there is here played by 
the Barnes zeta function. For a recent analysis on the dimensional dependence of the Casimir 
energy for scalar fields with Dirichlet boundary conditions and the electromagnetic field in the 
presence of a spherical shell see [^ , Q, where the space dimension D has been dealt with as a 
parameter, and results for (in principle) all values of real D have been obtained. 

The article is organized as follows. In the next section we briefiy recall the definition of 
Casimir energies in terms of zeta functions. In Sect. 3, we shortly describe the method and 
derive the formulas that are subsequently needed in the context of Casimir energy calculations 
1 26, 27]. In Sect. 4 we consider the case of a scalar field. For Dirichlet boundary conditions, 
the energy in dimensions D = 2 up to Z) = 9 is given there. The interior and the exterior 
regions are treated separately. Afterwards, the changes in the procedure needed for Robin 
boundary conditions are explained, and the corresponding formulas are derived. Given that 
the Casimir energy of the electromagnetic field is determined by using the Casimir energy of 
a scalar field satisfying Dirichlet boundary conditions (TE modes) and a scalar field satisfying 
Robin boundary conditions (TM modes), these forms constitute the basis for the electromagnetic 



case, and nothing else needs to be calculated, as will be later described in detail (Sect. 6). Before 
that, Sect. 5 is devoted to the spinor field. Local bag boundary conditions, as well as global 
spectral boundary conditions, are considered. In the concluding section, Sect. 7, a summary of 
our main results, as well as details on how our method is indeed able to yield arbitrary accurate 
results, are given. 



2 The Casimir Energy 

The Casimir energy of a quantum field $(t, x) inside a spherical shell is formally given by 



Eca 



1 



XI '^fc' 



(2.11 



(we set h = c= 1) with the one-particle energies lj^ = ^/X^ being obtained from 

- A(pk{x) = Xkcpkix), 



(2.2) 



fulfilling also suitable boundary conditions. The field operator is, in our case, A = d^ — A, 



and we have ^(t,x) 



-iult J.( ;^' 



x). The Laplacian A is the one defined inside or outside the 



D = {d + l)-dimensional ball B^ = {x S IR | ||x|| < R} and the fields (j){x) must satisfy 
appropriate boundary conditions at ||x|| = R. 

The Casimir energy as given by the formal expression ( |2.lD is ill defined and has to be 
regularised. In the ^-function regularisation procedure, one writes 



Ecas = \lJi^X{s-l/2)\s=Q 



(2.3) 
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Here, ^ is an arbitrary parameter with dimensions of mass to yield the correct dimension for 
all values of s, and C,{s) is the C-function corresponding to the operator A. In some cases, Ecas 
will be divergent and, as is known and will be seen later on, renormalisation ambiguities may 
remain. 

In order to calculate Ecas according to the previous definition, we need information on the 
zeta function C,{s) in a neighborhood of s = —1/2. As we are dealing with operators in flat 
space, but satisfying boundary conditions on a d-dimensional sphere (d = -D — 1, the boundary 
of the D-dimensional ball), the eigenvalues will be implicitly given as the zeros of a polynomial 
P{Z^,Zl) involving Bessel or Hankel functions, according to whether one is considering the 
internal or the external domain, respectively. We will denote the associated zeta functions by 
(^™*(s) and C^^^{s). The total Casimir energy will be the sum of the two terms, that is 



Eca 



1 



2s 



C^'is - 1/2) + C^'^s - 1/2) 



(2.5) 



With just a few modifications, which involve the phase of the zeta function (see |41] for precise 
details), all the considerations above can be extended to the Dirac operator. The basic construct 
turns out to be the zeta function of the square of the Dirac operator and one encounters a minus 
sign in Eq. (12? 



3 The method 



The method to be used here has been developed in the seminal papers |26, 21, ^] and permits 
to compute the C-function starting from the (indirect) knowledge of the eigenvalues through an 
implicit relation of the kind 

P(Z^, (iOnlR), K, i^nlR)) = , Xnl = ^ll , (3.1) 

where n,l > are the principal and azimuthal quantum numbers respectively. The degeneracy 
d{l) of the eigenvalues and the index ui of the Bessel functions depend on I and on the dimension. 
Their explicit forms are strictly related to the fields and the boundary conditions. 
The (^-function can be expressed as an integral in the complex plane, that is 



C(^) = E |§ / ^"''^ lnk-'"''P{Z,,{k), Zl^ik)) dk, Res>^, (3.2) 



where the open contour 7 has to be chosen to run counterclockwise and to enclose all strictly 
positive solutions of Eq. ( |3.lD . The additional factor k~^'^^ has been inserted in order to cancel 
the pole at the origin, which is important when deforming the contour in the next step. In 
this way 7 can also include the origin. Here 6 is a number which depends on the asymptotic 
behaviour of P at the origin: in our cases it will be ±1 for scalars, but for spin 1/2 with mixed 
boundary conditions it turns out to be ±2. 

For explicit calculation, it is convenient to write Eq. ( p.2[ ) as an integral on the real axis. This 
can be done by deforming the contour 7 to the imaginary axis and by making the substitution 
k -^ iy. In general one has to be careful when deforming the contour in that no poles in the 
plane Re A; > are hit. In fact, for Robin boundary conditions one has 

P{Z^, Zl) = aZ^{k) + kZlik) = {a- u)Z^{k) + kZ^^i{k) = , (3.3) 

which may have solutions for A; ^ IR too, if a > u. To avoid these cases, in the following we will 
consider a < uq only, uq corresponding to the smallest eigenvalue. 

With this assumption we can write the (^-function in the more simple form 



^ s_^^^ g ^(^) r y-2sd in[y-''-.p(^,, y)] dy , (3.4) 

71" t:^ Jo oy 



which is valid for 1/2 < Re s < 1 (for details see [^]). Here, P{v,y) = P{Zi,{y)^Zl{y)) is 
a polynomial like P (aside, possibly, from an irrelevant sign) and the Zy[y) = Zy{iy) are the 
modified Bessel functions corresponding to Z. In order to compute the Casimir energy we need 
the C" function at s = —1/2 and so we have to do an analytic continuation of Eq. (p.4|). 

With this aim, let us now employ the asymptotic expansion of the modified Bessel functions. 
For large values of i/, we have 
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where rj = \/l + z^ + ln[z/(l + \/T+z^)]. The first few coefficients Uk and Vk are hsted in [^, 
while higher order coefficients are immediate to obtain by using the recursion relations 

Uk+l{t) = h\l-t^)uUt) + lJ {l-5T^)Uk{T)dT, (3.9) 

Vk+i{t) = Uk+i{t) - -t{l - t^)uk{t) - t\l - t^)u'kit) , (3.10) 

1 yr^ , , 

As we shall see explicitly in the following, the above behaviour of Bessel functions permits 
us to write 

lnP{u,zu)r.\nFiu,z) + J2 ^, (3-12) 

71=1 

an expression which is valid for large values of v. The function F is related to the exponential 
factors in Eqs. (|3.5| )- (l3.8| ), while the coefficients Dn{t) are related to S^ and are polynomials in 
t. More precisely 

2n 
Dn{t) = J2 ^nfct"+'^ • (3.13) 

fc=0 

Note that when b = ±1, all Xnk with odd k vanish. Of course, -F, Dn and Xnk depend on the 
specific problem under consideration. We will specialize them for every case. 

The trick consists now in subtracting the asymptotic behaviour from the integrand function 
and in integrating the asymptotic part, with arbitrary s, exactly. We thus get 

N 

Cis) = Zois) + Z{s) + J2 Ms)- (3.14) 

n=-l 

Here, 

Zo(s) = 5(Z) - 2) (i[0] ^^^^^ Hdz z-^'-^ In P{0,z) (3.15) 

vr Jo oz 

is the contribution due to f/ = 0, which is present only in two-dimensions and has to be treated 
specifically for any case. Z(s) represents all the other terms with the asymptotic contributions 
subtracted, that is 



L P{ui,zui)- In F{ui,z)-J2^ 

[ n=l ^l 



dz 
and An are the integrals of the asymptotic part. They read 12' 



(3.16) 



Ms) = -Y{s)^m{s + n/2) J2 xnk^—^^ , n > 1 , (3.17) 

A_, +Ao = ^^^^ Y. ^(0 r ^^ ('""^y'^i- H{z^i)-'"''n^i, ^)] 

vr ^^^ Jo dz 

= c.i{s)Cm{s-1/2)+co{s)Cu{s), (3.18) 



Cxis) = E d{lK'' ■ (3.19) 



Eq. ( 3.16 ) is convergent for {D — 2 — N)/2 < Re s < 1, thus for our aim it is sufficient to 
subtract N = D asymptotic terms. This means that with N = D we can directly put s = —1/2 
in Eq. ( p.l6| ) and perform the integral numerically. 

As we shall see in the explicit examples, the base (^-function, C_yv', can be conveniently ex- 
pressed in terms of the Barnes zeta function |E3], defined as 



CB{s,a;d) = J2 7 — ] ] ] ^ = Zl ^"(^)(" + '^) ' ' 

^^ {a + mi + ...+mdr 



e„(d) 



m=o ^ ' ^ ""^ n=0 

(d + n-1)! 



n!(d-l)! ' 
for Re s > d. Obviously, there is an expansion of the kind 

d-l 

en{d) = J2gUd){a + nr 

Q = 

and this yields the expansion of the Barnes zeta function in terms of Hurwitz zeta functions 

mm 



d-l 

CB{s,a;d) = ^ga(d)CH{s-a,a). (3.20) 

For example, for d = 2, we trivially get 

CB{s,a;2) = Ch{s- I, a) + (1 -a)C//(s,a). 

One can show that the Qa {d) are connected with the generalised Bernoulli polynomials [^] . This 
allows to determine, in a direct way, the residues and finite parts of the Barnes zeta function of 



the problem at hand. As a result, the asymptotic contributions in ( 3.14 ) are readily computed. 

4 The scalar field 

The field equation for this case reads 

- A0fc(x) = Afc0fc(x), (4.1) 

and has to be supplemented with Dirichlet or Robin boundary conditions. Here, A is the 
Laplace operator inside or outside the D = (d + l)-dimensional ball and we impose Dirichlet 

(^(af)||^l=^ = 0) or Robin ([a(/)(x)| + (j)' {x)]\^^^=ji = 0) boundary conditions. 
In polar coordinates the solutions are 

with vi = I + [D — 2)/2, the fu{r) being Bessel functions and the ^z+z)/2(^) hyper spherical 
harmonics 



4.1 Scalar field with Dirichlet boundary conditions inside a spherical shell 

In this case, the fi, are Bessel functions of the first kind and thus the eigenvalues X^n = ^fn *^^ 
defined through 

and have degeneracies given by d{l) = [21 + d — 1) m/ ,~,-.| • From the last equation, it easily 



follows that |4|, I? I 



CM = CB[2s,-^;d)+CB(2s,^;d] . 



In this case, 

and, as a consequence, 

Fiu,z) 



P{v,z) = h{z), 
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(4.2) 



(4.3) 



(4.4) 
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The asymptotic contributions have been calculated to be 

1 r(.-i) 






4^ r(s + i) 



C^(s-l/2), 



1 



(,m{s). 



(4.5) 



(4.6) 



The C-function for the present situation is obtained by means of Eqs. ( 3.13 )-( p.l7 ) with the 
definitions above. 

As already anticipated in the previous section, in two dimension we have an aditional con- 
tribution that has to be computed explicitly. With this aim, we recall that, for large z, 



«^'^i5{'+^+°^" 



and thus we can write 
^o(s) 



smlvrs 



vr 



{|</..-|ln/„M 
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1 

'ii"^ 2(s-l/2) 16(s + l/2) 
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(4.7) 



where the poles at s = ±1/2 are shown explicitly. The integrals are now convergent for s = —1/2 
and can be computed numerically. 



4.2 Scalar field with Dirichlet boundary conditions outside a spherical shell 



Now the radial part of the solutions are Bessel functions of the third kind (Hankel functions), 
while f/ and d{l) remain the same. Thus, we have 



ui = l + 



D-2 



(4.8) 



d{l) = {2l + d-l) 



{l + d-2)l 

u{d-iy. 



Cat = Cb ('^s. 



d-l 



^±l-d)+CB(2s,^^-d 



(4.9) 
(4.10) 



P{u,z) = K,{z), 



F{u,z) 



-UTj 



1 + ^2)4 



(4.11) 
(4.12) 



InS, 
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(4.13) 
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Ca^ {s - 1/2) , 



(4.14) 
(4.15) 



Owing to the particular relation between Si and S3, the coefficients Dn{t) differ from the 
corresponding coefficients one has in the internal case just for the trivial factor (—1)"'. The 
same holds also for the quantities ^^(•s). 

In two dimensions we have to consider also the contribution due to z^ = 0, which can be 
obtained with the same arguments as in the previous case, Eq. ( |4.7l ). The result is 



Zo{s) 



sm vrs 



/oo 
dz z 
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1 f} 

dz z-^'—lnKoiz) 
dz 



l^"^°(^) + ^ + ^-8i^ 
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1 



4s 2(s-l/2) 16(s + l/2)J' 



(4.16) 



The numerical results corresponding to the ^-functions inside and outside the shell and the total 
Casimir energy are reported in Table 1 for the choices D = 2,. ..,9. For the interior space for 



D = 2 and Z? = 3 as well as for D = 3 and the exterior space, our results agree with |22|. For 
the whole space in D = 3 the result is given in |23, 3£]. 



D 


C(-l/2) Inside 


C(-l/2) Outside 


Casimir Energy 


2 


+0.0098540 - 0.0039062/e 


-0.0084955 - 0.0039062/e 


+0.0006793 - 0.0039062/e 


3 


+0.0088920 + 0.0010105/e 


-0.0032585 - 0.0010105/e 


+0.0028168 


4 


-0.0017939 + 0.0002670/e 


+0.0004544 + 0.0002670/e 


-0.0006698 + 0.0002670/e 


5 


-0.0009450 - 0.0001343/e 


+0.0003739 + 0.0001343/e 


-0.0002856 


6 


+0.0002699 - 0.0000335/e 


-0.0000611 - 0.0000335/e 


+0.0001044 - 0.0000335/e 


7 


+0.0001371 + 0.0000214/e 


-0.0000555 - 0.0000214/e 


+0.0000408 


8 


-0.0000457 + 5.228 x lO'Ve 


+0.0000101 + 5.228 x IQ-Ve 


-0.0000178 + 5.228 x IQ-^/e 


9 


-0.0000230 - 3.769 x lO^V^ 


+0.0000094 + 3.769 x IQ-Ve 


-0.0000068 



Table 1: Scalar field with Diriciilet boundary conditions. Values of the zeta function at s = 1/2, 
for the inside and the outside regions of a spherical shell, and values of the Casimir energy. The presence 
of the cutoff e for all even dimensions is to be noted. In such cases, the Casimir energy is divergent and 
needs to be renormalised. 



4.3 Scalar field with Robin boundary conditions inside a spherical shell 

In the case of Robin boundary conditions the radial part of the solution is a combination of 
Bessel functions with derivatives. For the interior case we have Bessel functions of the first kind 
and their eigenvalues are determined through 



f 1 - y -/3j Jui{^l,n) + ^l,nJlii^l,n) = 0. 



(4.17) 



Here we have put a = 1 — D/2 — (3 and, in the spirit of Sect. 0, we have to restrict ourselves to 
the case /3 > 1 — D/2 — uq. The choice /3 = represents Neumann boundary conditions. 
Also for this case vi, d{l) and Ca^ are given by Eqs. (|4.8|)-(4!Tc|), now with 



P(u,z)=(l-^-P]l,iz) + zlUz), 



F{v,z) 



^e^vil + z^)\ 
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1 r(.-i) 

A^ r(s + i) 
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In two dimensions we have to consider also the contribution 
Zq{s 

TV UO 
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D 


C(-l/2) Inside 


C(-l/2) Outside 


Casiniir Energy 


2 


-0.3446767 - 0.0195312/e 


-0.0215672 - 0.0195312/e 


-0.1831220 - 0.0195312/e 


3 


-0.4597174 - 0.0353678/e 


+0.0120743 + 0.0353678/e 


-0.2238215 


4 


-0.5153790 - 0.0447159/e 


-0.0060394 - 0.0447159/e 


-0.2607092 - 0.0447159/e 


5 


-0.5552071 - 0.0489213/e 


+0.0030479 + 0.0489213/e 


-0.2760796 


6 


-0.5949395 - 0.0513727/e 


-0.0128321 - 0.0513727/e 


-0.3038858 - 0.0513727/e 



Table 2: Scalar field with Neumann boundary conditions (or Robin with the choice /? = 0). 

Values of the zeta function at s = —1/2, for the inside and the outside regions of a spherical shell, and 
corresponding values of the Casimir energy. 



4.4 Scalar field with Robin boundary conditions outside the spherical shell 

As for Dirichlet, the only difference between the interior and the exterior case consists in the 
replacement of Bessel functions with Hankel functions. Eqs. ( |4.8| )-( pITo| ) are valid again, while 



P{u, z) = {l--^-l3^ K,{z) + zKl{z) , 
F{v,z) = ^e-^^{l + z^)\, 



ln,1^^^..3-.0-E^ 
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A_,(») 
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Cm {s - 1/2) , 



Ms) = -^Cm{s). 

For the v = Q contribution, we have in this case 
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ln{aKo{z) + zK^{z)) 



— HaKoiz) + zK^{z)) + 1 - 1 + g - a") 1 



1 



4s 2(s - 1/2) 



a 
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All numerical results corresponding to Neumann boundary conditions (or Robin ones with (3 = 0) 
are exhibited in Table 2. For D = 2 the result is given in p2|, for Z) = 3 in [23|. 



5 Spinor field on the D-dimensional ball: bag boundary condi- 
tions 



We now consider spinor fields, see [48, 37]. The eigenvalue Dirac equation on the Euclidean 
D-hal\ is 



iF^V^Vi = ±A:V±, F^'^F'^) = <7'^^ 



(5.1) 



10 



and the nonzero modes are separated in polar coordinates, ds^ = dr'^ + r'^dO?^ in standard 
fashion to be regular at the origin (C and A are radial normalisation factors), 






r 









Here the Z^ (O) are well-known spinor modes on the unit [D — l)~sphere (some modern refer- 
ences are [^, |5^, ^) satisfying the intrinsic equation 

- i7' V,zi") = ±A„zi"\ (5.3) 

where 

D-l 

A„ = nH — , n = 0, 1,... . 

For D > 2, each eigenvalue is greater than or equal to 1/2 and has degeneracy 

The dimension, dg, of V'-spinor space is 2^/^ if D is even. For odd D it is 2^^+^)/^ and has been 
doubled in order to implement the boundary conditions. The projected 7-matrices are given by 

5.1 Spinor field inside a spherical shell: bag boundary conditions 

For bag — also called mixed — boundary conditions, we apply P+iJj = at r = 1, where the 
projection is given by 

P+ = i(l-ir5r^n^), (5.5) 

in terms of the inward normal n^. For the geometry of the ball 

1/1 a~ 



and so for ipj. , 



and for ip^ , 



^+ 2V-^1 1 



Jn+D/2{k) — TJu+D /2-l{k), 



Jn+D/2-lik) = TJn+D/2ik), n = 0, 1, 2, . . . . 

Thus, taking f„ = n + (D — 2)/2, the implicit eigenvalue equation is as in [|^ 

J^ik)-J^+,ik) = 0, (5.6) 

while the degeneracies are 

din) = ds(^^'^-^). (5.7) 

11 



In two dimensions the degeneracy is just 2. 

In summary, all the relevant functions for this case are 



f„ = n + 



D-2 



d{n) = d. 



{n + D-2)\ 
n\{D-2)\ ' 



d{n) = 2 , for D = 2 



(5.8) 



Q^f{s) = dsQB{2s,D/2-l■d) 



Cm{s) = 2Cr{2s) , for D = 2 



P{v,z) = ll{z)+Il^,{z), 



F{v,z) 



(I_t)e2^'?(l + z2)l 



■KUZ'^ 
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1 



2(1 -t) 



S? + Si-2iSiS2 l-E^ 



n=l 



A-i{s) 



^o(s) 



1 r(.-i) 

2^^ r(s + i) 
1 r(s + i 



2^^ r(s + i) 



C^(s-l/2), 
CAr(s), 



with ^^ the Riemann ^-function. 

The contribution of i^ = 0, which we have in two dimensions reads here 



Zo{s) 



""("^)'[°l^^'d..--#ln(/o^(.) + /f(z)) 



/oo 
1 



2s 



dz 



1 



+ 



1 



2s (s-1/2) 8(s + l/2) 

5.2 Spinor field outside a spherical shell: bag boundary conditions 

As in the scalar cases, we must simply replace Bessel with Hankel functions. Eqs. (|5^ 
( |5.9| ) provide some quantities needed in the computation, while for the rest, we get 

P{v,z)=Kl{z)+Kl^^{z), 



and 



F{v,z) 



4z^(l + t)e-2'^''(l + z2)| 



■KZ^ 



In 



^ ^S^ + Sl + 2tS3S4l)~f:^ 



2(1 + i) 



n=l 
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D 


C(-l/2) Inside 


C(-l/2) Outside 


Casiniir Energy 


2 


-0.0058312 + 0.0078125/e 


+0.0213677 + 0.0078125/e 


-0.0077683 - 0.0078125/e 


3 


-0.0605944 - 0.0050525/e 


+0.0198217 + 0.0050525/e 


+0.0203863 


4 


+0.0059074 - 0.0028381/e 


-0.0101965 - 0.0028381/e 


+0.0021445 + 0.0028381/e 


5 


+0.0250447 + 0.0025110/e 


-0.0089912 - 0.0025110/e 


-0.0080268 


6 


-0.0030244 + 0.0011715/e 


+0.0046183 + 0.0011715/e 


-0.0007969 - 0.0011715/e 


7 


-0.0108618 - 0.0011745/e 


+0.0040247 + 0.0011745/e 


+0.0034186 



Table 3: Massless spinor field with mixed boundary conditions. Values of the zeta function at 
s — —1/2, for the inside and the outside regions of a spherical shell, and values of the Casiniir energy. 



A_,{s) 



Ao{s) 



1 



r(,-i) 



2^ r(s + i) 
1 r(. + i) 



Cat {s - 1/2) , 



20F r(s + l 
In the same way, for z^ = we obtain 

sin(7rs)d[0] '' '"i 



Cm{s). 



Zo{s) 



/oo 
dz z 

1 



dzz--^^^HKl{z) + Kl{z)) 
oz 

d ■ '^^2,.^ , r.2^ ,, _ 1 1 



^Jn{K^{z)+Kt{z))+2 + - + ^^^ 
1 1 



2s (s-1/2) 8(s + l/2) 

The numerical results for spin 1/2 with bag boundary conditions are given in Table 3. The 
D = 3 result is the one found already by Milton [15| (albeit with far less precision). 



5.3 Spinor field with global spectral boundary conditions 



We shall now obtain the results for spectral boundary conditions [[48| , 37|. Such boundary 
conditions are imposed by setting equal to zero, at r = 1, the negative (resp. positive) Z-modes 
of the positive (resp. negative) chirality parts of ip, the rest of the modes remaining free. 

Roughly speaking, spectral conditions amount to requiring that zero-modes of (|5.1| ) should 
be square-integrable on the elongated manifold obtained from the ball by extending the narrow 
collar (of the approximate product metric dr'^ + dJ7^) just inside the surface, to values of r 
ranging from 1 to oo. This will be so if the modes of A = FT^Va with negative eigenvalues are 
suppressed at the boundary (e.g. [5^] -[pi]]). 

From (5^) and (|5.3D , the boundary operator is Aq = F'"F"Va|r=i ^^'^ ^^^ ^'^^ eigenstates 



^n 



Z 



An) 



An 



(n) 

+ 
(n) 



An 



(n) 
(n) 



(n) 
(n) 



(5.9) 



Thus, from ( |5.2D we see that the negative modes of Aq are associated with the radial factor 
Jn+D/2-i{^f)- Taking u as before, v = n + {D — 2)/2, the implicit eigenvalue equation reads 

Ju{k) = 0. 
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D 


C(-l/2) Inside 


C(-l/2) Outside 


Casiniir Energy 


2 


-0.0093152 + 0.0319762/e 


+0.0100172 + 0.0319762/e 


-0.0003510 - 0.0319762/e 


3 


-0.1710212 -0.0037705/e 


+0.0019763 + 0.0037705/e 


+0.0845225 


4 


+0.0082635 - 0.0118316/e 


-0.0040473 -0.0118316/e 


-0.0021081 + 0.0118316/e 


5 


+0.0680217 + 0.0019471/e 


-0.0009007 - 0.0019471/e 


+0.0335605 


6 


-0.0042224 + 0.0049069/e 


+0.0017603 + 0.0049069/e 


+0.0012311 - 0.0049069/e 


7 


-0.0290717 - 0.0009256/e 


+0.0003983 + 0.0009256/e 


+0.0143367 


8 


+0.0020298 - 0.0021417/e 


-0.0007907 - 0.0021417/e 


-0.0006196 + 0.0021417/e 


9 


+0.0128994 + 0.0004353/e 


-0.0001787 - 0.0004353/e 


-0.0063604 



Table 4: Massless spinor field with global spectral boundary conditions. Values of the zeta 
function at s = —1/2, for tlie inside and the outside regions of a splierical shell, and corresponding values 
of the Casimir energy. 



The degeneracy for each eigenvalue is 
d{n) = 2ds 



n + D-2 
D-2 



d{n) = 4 , for D = 2 . 



The relevant boundary zeta function reads now 

CAs) = 2dsCB{2s, D/2 -l;d), Ca/-(s) = 4Cij(2s) , for D = 2 



(5.10) 



(5.11) 



As we see, apart from the degeneracy of the eigenvalues and the relation between (!^_\f and the 
Barnes (^-function, the rest of the argumentation is identical to that for the scalar case with 
Dirichlet boundary conditions. Thus, Eqs. ( [4.3D -(p^) remain valid once the above definitions 
are used. 

For the exterior space we have to employ Eqs. ( ^.10 ) and (5.11) in Eqs. ( [4.1lD -( ^.16| ), but 
it has to be noted that here, in contrast with the interior case, I'l = I + D/2, as a result of the 
normal vector changing its sign. This means that there is no z^; = contribution. 

The numerical results for this case are listed in Table 4, for D = 2, ...,9. 



6 Electromagnetic field in a perfectly conducting spherical shell 

The Casimir energy of the electromagnetic field is, essentially, the sum of a Dirichlet and of a 
Robin scalar field (with a specific value for /3, see Eq. ( 4.17| )), the only difference being that the 
angular momentum / = is to be omitted. An exception is -D = 2, where the vector Casimir 
effect consists of only the transverse magnetic mode contributions. Being precise, in the interior 
of the shell one has for the transverse electric (TE) — respectively for the transverse magnetic 
(TM) modes — the following boundary conditions [ffl, |62|, 



{D/2 - l)Ji,^{uji^nr) + ^l,nJ'ui{^l,nr) 



r=R 



for TE-modes, 
: 0, for TM-modes. 



(6.1) 



The condition for the TM-modes is of Robin type with j3 = 2 — D. Since the I = mode has 
to be omitted, the minimum eigenvalue in this case is /xi = D/2 and therefore we can apply 
the method for any [5 = 2 — D>1 — D. Thus, in order to get the Casimir energy of the 
electromagnetic field, we must simply repeat the computation of Sect. ^ for Dirichlet and Robin 
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D 


C(-l/2) Inside 


C(-l/2) Outside 


Casimir Energy 


2 


-0.3446767 - 0.0195312/e 


-0.0215672 - 0.0195312/e 


-0.1831220 - 0.0195312/e 


3 


+0.1678471+0.0080841/6 


-0.0754938 - 0.0080841/e 


+0.0461767 


4 


+0.5008593 + 0.0231719/e 


-0.1942082 -0.0564056/e 


+0.1533255 - 0.0332337/e 


5 


+1.0463255 + 0.1838665/e 


-0.2981425 - 0.1838665/e 


+0.3740915 



Table 5: Electromagnetic field in a perfectly conducting spherical shell. Values of the zeta 
function at s = —1/2, for the inside and the outside regions of a spherical shell, and corresponding values 
of the Casimir energy. Note that in even dimensions, in contrast with the scalar field, the divergences 
arising from the inside and the outside energies are different. This is due to the fact that (only in even 
dimensions) the / — mode explicitly contributes to the poles of the ^-function, such contribution being 
absent in the scalar case. 



boundary conditions with (3 = 2 — D and add them up. We have to exclude everywhere the 
/ = mode and this means that also the base C-function is slightly modified, in the way 

-2s 



C^(s) = Ce ( 25, ^; dV Cb ("2./ 



1 



d 



d-l 



(6.2) 



The results for the electromagnetic field are summarised in Table 5. D = 2 \s the Neumann 
result, D = 2> is the well known figure first obtained by Boyer |li and later recalculated in Refs. 



7 Discussion and conclusions 



In Ref. |26|| , two of the present authors developed a new, seminal approach for finding represen- 
tations of the zeta function associated with the Laplace operator on the D-dimensional ball. At 
that stage, dimension by dimension was considered, but soon a refined and generalised technique 
was provided in subsequent works [^, ^ . Making use of the Barnes zeta function |^^ , dimen- 
sion can easily be dealt with as a parameter and several different fields can also be treated on 
the same footing. The representations derived are valid for all values of the complex parameter 
s and it is up to the practitioner's needs or wishes at which values of s the zeta function is to be 
evaluated. In previous work our concern was of a more mathematical nature and we considered 



function values and residues appropriate to find heat-kernel coefficients |26, 27, 37], as well as 
the derivative at s = |3C]. Since then, a number of proclaimed "new" methods have been 
developed in the literature. 

Our aim in the present article has been to show explicitly that Casimir energies for the 
big family of the more usual configurations can be obtained in fact from general formulas, also 
in quite non-trivial situations, where the boundaries are not flat plates, the flelds are spinorial 
(rather than scalar), and also when the boundary conditions are very general and rather involved. 
We have gone much beyond previous work in that we are not here restricted any more to a very 
specific field in a specific dimension with a specific boundary condition, but give general formulas 
for basically any possible situation that can arise in practice, involving spherically symmetric 
boundaries. 

Some comments on the precision and accuracy of the numerical procedure employed are in 
order. It is clear from the analysis developed in the previous sections, that a numerical evaluation 
of the asymptotic terms Ai{s) to any desired accuracy is immediate, using the formulas given 
there. These contributions are always represented by known special functions and using available 
programs, such as Mathematica, the accuracy with which these are calculated is readily obtained. 
Imposing accuracies of, e.g., 10"^*^ or more, we get results in negligible cpu time. 
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The only problem (if at all) with the numerical analysis is the computation of Z{s), Eq. ( |3.14| ). 
It is twofold. On one hand, the integration, up to infinity, of the combination of Bessel functions 
is not strictly possible, using the exact form of the Bessel functions. On the other hand, the 
angular summation, up to infinity, cannot be performed exactly. For large angular momenta, the 
Bessel functions take a rather complicated form, what renders exact summation not possible. For 
that reason, the following procedure has been applied throughout (the error bounds given below 
are for Dirichlet boundary conditions, but very similar relations hold for the other conditions 
considered). 

We have dealt with the infinite integration as follows. The main contributions always come 
from small values of z, and thus we split the integral into 

B roo 

dz + dz. 
Jb 

Whereas in the first integral the Bessel functions themselves are used for the integration, in the 
second integral their asymptotic expansion for large arguments is employed. The value of B is 
computed with the help of an adaptative procedure, such that the integrand and its asymptotic 
expansion differ, at B, by less than, say 10~^^. Typically i? = 10 is already sufficient. Given that 
the asymptotics of the Bessel functions is a simple polynomial in powers of (1/z), the integration 
up to infinity is very easily done. 

Let us now assume that the contribution of the first L angular momenta has been calculated 
as described. In order to get a numerical approximation for the angular momentum sum, from 
L + 1 to infinity, we proceed as follows. The idea is that, for sufficiently large values of L, 
the integrand can be replaced by its uniform asymptotic expression. For Dirichlet boundary 
conditions this amounts to going through the following steps: 

"' Dn{t) 



1 . . /"CO 

/ poo \ °° 

/ dzDN+iit)) J2 dil) 

^^0 ^ l=L+l 

/ roo \ °° 

+ ^ dzDN+2{t)j E ^W 



In I,{uz)- In ^=^ --Y^^ 

"^ ' V2^(l + z2)l/4 ^^ yn 



-N 



u~^-' + 



l=L+l 



+ (^1^" dz Div+2(t)) (cm{{N + l)/2) - Y. d{l),.-^-A + .. 



(7.1) 



Again, the integrals over the uniform asymptotics are simple and can be performed analytically. 
In this way, a closed expression for the approximation is found, the value of L being again 
determined by an adaptative procedure. By definition, the difference Z£"* — ZY^K equals the 
contribution coming from I = L. The value of L is determined such that the difference Zf^^ — 
ZY+n obtained from ([7.1|), agrees up to say lO"^'' with the contribution from I = L calculated 
previously. Depending on the dimension, the values of L range from 6 (for D=9) to 49 (for 
D=3). 

In summary, as explained, this procedure takes fully into account the integrals of infinite 
range as well as the summation up to infinity. The error bounds can thus be imposed at will in 
the single steps and this guarantees that the results given are always numerically precise, up to 
any pre-established digit. To our knowledge, this does not apply to any other method. 
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In the cases when partial results were known, we have compared our numbers with those 
while improving always, by several digits, such known values and deriving, for the first time, 
a lot of new ones, for different fields (e.g. results for the exterior space in the case of the 
electromagnetic field) and different boundary conditions (e.g. for spectral boundary conditions, 
and for bag boundary conditions in any dimension). For the scalar field with Dirichlet boundary 
conditions we have re-obtained, in particular, the known result that for D even the energy is 



divergent |39|. Here it still remains unclear if there may be a natural way to get, unambiguously, 
a finite answer with physical sense for this case. In odd dimension, D = 2n — 1, the sign of the 
Casimir energy seems to be determined by the sign of (—1)". For even dimension, D = 2n, one 
also finds the alternating structure (—1)""*"^ for the finite part of the Casimir energy; however, 
its interpretation is unclear due to the presence of the pole. Similar comments hold for the 
interior and exterior contributions separately, with the same problems of interpretation. For 
Neumann boundary conditions, in all the dimensions calculated, the Casimir energy is negative. 
Similarly, one can describe the results summarized in Tables 3-5. In all cases we have been able 
to obtain general, highly accurate expressions which, by fixing some parameters, provide us with 
the desired specific example and yield a numerical answer of arbitrary precision (just by adding 
the convenient number of terms of the corresponding series). 

Disappointing as the mentioned — quite well known — ambiguities may be (specialists in the 
field are quite used to them by now), even more so is the fact that no general pattern seems 
to arise from our general formulas which might hint towards the physical understanding of the 
final sign of the energy. Here we have been able to demonstrate, without reasonable doubt, the 
existence of the two classes of Casimir force: attractive and repulsive, but are unable to give the 
rule according to which one or the other will show up at a particular instance. Further insight 
will be needed to clarify this point. 
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